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The Sombrero galaxy (NGC 4594) is an Sa galaxy with a symmetric dust ring. We have used the Large APEX BOlometer CAmera 
(LABOCA) at 870 i^im and the MAx-Planck Millimeter BOlometer (MAMBO-2) at 1.2 mm to detect the dust ring for the first time at 
submillimetre and millimetre wavelengths. We have constructed a model of the galaxy to separate the active galactic nucleus (AGN) 
and dust ring components. The ring radius at both 870 f^im and 1.2 mm agrees well with the radius determined from optical absorption 
and atomic gas studies. The spectral energy distribution of the ring is well fitted by a single grey-body with dust emissivity index 
/? = 2 and a dust temperature r;j = 18.4 K. The dust mass of the ring is found to be 1.6 ± 0.2 x 10^ Mq which, for a Galactic gas-to-dust 
ratio, implies a gas mass that is consistent with measurements from the literature. 
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> 1. Introduction 



The Sombrero galaxy is one of the best studied early-type spi- 
ral galaxies. It is an Sa galaxy with a prominent dust lane, 
which is actually a symmetric dust ring seen nearly edge-on, 
and a nucleus that h osts a lO^M© supermassive black hole 
dKormendv et al.l [1996). These two features are the highest 
surface brightness infrared sources in the galaxy (iBendo et al.l 
l2006ai hereafter B06). At a distance o f 9. 4 Mpc (the averag e 
of measurements from iFord et al.|[T996l and iTonrvetai] 1200 ih . 
it is an ideal laboratory for investigating different astrophysical 
processes, in particular, the connection between star formation 
and an active galactic nucleus (AGN) and the effects of these 
mechanisms on nuclear and global galaxy properties. The galaxy 
is also ideal for studying the separate spectral energy distribu- 
tions (SEDs) of diffuse interstellar dust and the dust heated by 
an AGN, since the majority of diffuse dust is located in a ring 
relatively far from the AGN, making it thus relatively simple to 
model either component (B06). 

Until recently, investigation of this gal axy in the far- 
infrared (FIR ) and submillimetre (submm) (e.g. lRice et al.ll998t 
iKrause et all 12006) was limited to the global SED due to the 
resolution of IRAS. B06 investigated the mid- and far-infrared 
SED of the separate components using Spitzer (3.6-160 i^im) and 
SCUBA (850 nm) observations. However, only the AGN compo- 
nent was detected with SCUBA, and the SCUBA map covered 
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only ~ 3' X 2', while the diameter of the dust ring as seen by 
Spitzer measures ~5' (B06). 

We present observations of the Sombrero galaxy 
(NGC 4594) at 870 ti m with the Larg e APEX BOlometer 
CAmera (LABOCA; Krevsa et all l2003t ISiringo et all l2007h 
and at 1.2 m m with the MAx- Planck Millimeter BOlometer 
(MAMBO-2: lKrevsa et al.ll998l) . In order to separately examine 
the dust ring and the AGN, we have modelled the dust emission 
using the method described in B06. We present the first results 
from this analysis, and focus on the FlR/submm emission from 
the dust ring. Investigation of the FlR/submm emission from the 
nucleus will be the subject of a separate paper (Vlahakis et al. in 
prep.). 



2. Observations and Data Reduction 

The 870 ^im data were taken with the newly commissioned 
LABOCA bolometer array, located on the Atacama Pathfinder 
Experiment (APEX) telescope on Chajnantor, Chile, during its 
science verification phase in May 2007. The full width half max- 
imum (FWHM) of the point spread function (PSF) at 870 (,im is 
~19". The obser vations were carrie d out using a raster-spiral 
observing mode dWeiss et al.l l2007l) in which each scan con- 
sists of 4 spiral pointings offset by 27", each spiral providing 
a fully sampled map of the full field of view (FOV) of LABOCA 
(11.4'). The total map thus covers an area slightly larger than 
the FOV. The total on-source integration time was 40 minutes, 
reaching an rms level of ~8 mJy beam"' in the unsmoothed map. 
We reduced the data with the BoAQ software package. The data 
reduction process included correction for atmospheric zenith 
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Fig. 1. Submm/mm images of the Sombrero galaxy. Left panel: 870 ^im LABOCA map, smoothed to a resolution of ~22", shown 
as a contoured colour scale (top) and as 870 (,im contours overlaid onto an optical image (bottom). The rms in the smoothed map is 
3.4 mJy beam"^ Contours show flux densities of 6.8, 10.2, 17, 34, 85, and 136 mJy beam"'. The central 8' x 5' of the map is shown. 
Right panel: 1.2 mm MAMBO-2 map, smoothed to a resolution of ~15", shown as a contoured colour scale (top) and as 1.2 mm 
contours overlaid onto an optical image (bottom). The rms in the smoothed map is 1.1 mJybeam"^ Contours show flux densities of 
2.2, 3.3, 11, 33, and 66 mJybeam The central 8' x 5' of the map is shown; we note that the peaks at the edges of the image are 
noise features towards the edges of the map. 



opacity (opacity, t, values ranged from 0.2 to 0.3), flat-fielding, 
despiking, coiTelated sky noise removal, and the removal of fur- 
ther correlated noise due to instrumental effects. The data were 
calibrated using observations of Mars and Uranus. The uncer- 
tainty due to calibration is the main source of uncertainty, and is 
estimated to be ~IQ%. For presentation purposes, we smoothed 
the final map to a resolution of ~22" by convolving it with a 
Gaussian. 



We carried out the 1.2 mm observations with the MAMBO- 
2 117-element bolometer array on the 30m telescope on Pico 
Veleta, Spain, during the pooled observing sessions in December 
2006 and March 2007. The FWHM of the PSF at 1.2 mm is 
~11". We used a standard on-the-fly mapping technique and 
a wobbler switching mode with a chop throw of 240" to map 
an area of ~ 9' x 5'. The total integration time was 96 min- 
utes, reaching an rms of ~3 mJy beam ' in the unsmoothed map. 
Opacity (t) values ranged from <0.1 to 0.3, while the sky noise 
was consistently low. We reduced the data using standard meth- 
ods in MOPSIC, part of the GILDAS software package. This in- 
cluded correction for atmospheric opacity, flat-fielding, despik- 
ing, and removing correlated sky noise. The data were calibrated 
using observations of standard calibration sources, resulting in a 
calibration uncertainty of ~15%. For presentation purposes, we 
smoothed the final map to a resolution of ~ 15" by convolving it 
with a Gaussian. 



3. Results 

The 870 nm LABOCA image and the 1.2 mm MAMBO-2 im- 
age are shown in Fig. [T] Contours showing levels of equal 
flux density are overlaid onto the images. These same contours 
are also shown overlaid onto an optical B-band image (from 
the Spitzer Infrared Nearby Galaxies Survey (SINGS) Legacy 
Data; Kennicuttet al. 2003). The dust ring and the nucleus are 
clearly detected in the LABOCA map. While the nucleus was 
detected previously at 850 \im with SCUBA and at 870 t-im with 
the Heinrich Hertz Telescope (B06; iKrause et al.ll2006l respec- 
tively), here the dust ring is detected for the first time at submm 
wavelengths. The dust ring and nucleus are also detected in the 
MAMBO-2 map. In both the submm and mm images, there is 
no evidence of dust outside the ring or outside the plane of the 
ring, and we do not detect any substructure in the ring at this 
signal-to-noise level. 

In order to separate the nuclear and ring components, we fit- 
ted a model to the unsmoothed data using the method described 
in B06. To summarise this method, the model consists of two 
components, a Gaussian ring and a central point source, which 
were convolved with a Gaussian PSF with a FWHM equivalent 
to that for LABOCA and MAMBO-2; the fit was performed over 
a region corresponding to the optical disk. The image models are 
shown in Fig.|2]together with the residuals from the model fits. 

The 870 [,im and 1.2 mm flux densities for the two compo- 
nents are given in Table [T] together with the global flux density. 
The uncertainties given are those resulting from the model fit- 
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Fig. 2. Image model (top) and residual image (bottom) made by subtracting the model image from the observed image for the 
LABOCA 870 (xm data (left panel) and the MAMBO-2 1 .2 mm data (right panel). 



Table 1. 

NGC 4594. 



Flux densities for the separate components of Table 2. Ring model parameters from fits to data. 



Wavelength 




Flux density (Jy) 


(urn) 


Global 


Nucleus/AGN Ring 


1200 


0.442 


0.185±0.002'' 0.257±0.006 


870 


0.924 


0.242±0.005 0.682+0.016 


870 




0.230+0.035* 


850 




0.25±0.0^ 



"Uncertainties are from the fits to the data and do not include cali- 
bration uncertainties (see Section|2j. 

*Previous submm measurements (B06: lKrause et al]|2006l) . 



ting process, while the largest source of uncertainly comes from 
the calibration. We note that the 870 [xm flux density we find for 
the nucleus is in excellent agreemen t with previous submm ob- 
servations (B06:^ Krause et al]|200 6'), which are listed in Table[T] 
for comparison. The ring model parameters are given in Table |2] 
along with the parameters determined from the 5.7-70 [xm data 
by B06. 

We have fitted a modified Planck function ('grey-body') to 
the 70|,im-1.2mm SED of the dust ring, using the 70- and 
160- [xm Spitzer measurements from B06 and our 870- \\m and 
1.2 mm measurements from Table [T] In this fit, the Planck func- 
tion was modified by a A'^ emissivity function with dust emis- 
sivity index fi-l, and the dust temperature was allowed to 
vary. The SED of the dust ring is found to be well fitted by a 
single grey-body with 7'd=18.4K (Fig.|3]l. 



Parameter" 


1.2 mm 


870 [.im 


5.7-70 urn^ 


Min/Maj axis ratio 


0.0757+0.0014 


0.0710+0.0013 


0.0990+0.0015 


Position angle 


90.43° ±0.04° 


90.69° +0.09° 


89.25°+0.06° 


Ring radius 


169.9"±1.0" 


162.0"+3.0" 


144.7"+0.5" 


Ring width 


33.0"±2.0" 


45.7"+0.6" 


22.9" + 1.1" 



"Definitions of the parameters are as given by B06. 
''Values from B06. 



We estimate the dust mass of the ring with 



(1) 



where ^gvo is the 870|,im flux density; D is the distance 
(9.4 Mpc); B(Td)9,n) is the Planck function for the fitted temper- 
ature T^; and k^-jq is the dust mass opacit y coefficient at 870 jxm. 
Assuming a value k^jq - 0.041 m^ kg"^ dLi & Drainell200ll) . we 
obtain Md - 1.6 + 0.2 x 10^ Mq. Here, the uncertainty is de- 
rived from the uncertainty on the 870 ^m flux density, which 
is the main source of uncertainty for the measured parameters. 
We note, however, that a potentially larger source of uncertainty 
could arise due to the value of k, which is notoriously uncertain. 

4. Discussion 

The dust ring detected at 870 (im and at 1.2 mm is wider than 
the dust ring detected at 5.7-70 ^im. This indicates that the dust 
ring is not monochromatic; the dust on the inside of the ring 
is warmer and radiates primarily at shorter wavelengths than the 
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Fig. 3. SED of the dust ring. A single grey-body with (3-2 is 
fitted to the 70 and 160 (im (Spitzer), 870 (xm (LABOCA), and 
1.2 mm (MAMBO-2) flux densities. The fitted dust temperature 
is 18.4K. 



dust on the outside. This interpretation is consistent with the dust 
temperature gradient sugge sted by the analysis of dust opacity 
measurements presented by lEmsellerd d 1 995.) . These results im- 
ply that the primary heating sources for the ring are located on 
the inside of the ring. 

The ~170" radius of the dust ring found at 870 ^m and 
1.2 mm is consistent with the outermost of t wo rings deter- 
mined in optical absorption (e.g.. lDettmajl986l : IWainscoat et alj 
Il990t lEm sellem 1995) and neutral hydrogen (Hi) emission 
dBajaia et al.l 11984 ). The 5.7-70 (im dust, molecular gas, ra- 
dio synchrotron emission, and part of the neutral hydrogen 
(Hi) emission lies in a -140" ring dBaiaia et alJll988L I199U 
iKrause et alj|2006!) . This ring is consistent with the innermost 
of the two rings detected in optical absorption. Thus, it appears 
that emission from wavebands more commonly associated with 
star formation originates from the ~140" ring whereas 'quies- 
cent' dust and gas emission, not necessarily associated with star 
formation, originates from the ~170" ring. 

The 1.2 mm data point is useful for constraining the 
Rayleigh- Jeans (RJ) tail of the SED, yet we note that the 70-, 
160-, and 870- (xm data alone are al so well fitted by this g rey- 
body. This result is in contrast with iGalliano et alj ( l2005h and 
iBendo et all (i2006b). who found that the RJ side of the SED de- 
viates from thermal dust emission at wavelengths longer than 
850 [im. 

The dust mass we calculated from our measured 870 ^m 
flux density and fitted dust temperature is about a factor of two 
higher than the d ust mass estimated from optical absorption by 
lEmseUeml d 1995b (8.3x10'' Mq, scaled to 9.4 Mpc). Part of this 
discrepancy can probably be explained if the re is an iriconsis- 
tency between the dust opacity model used bv^ Emselleml dl995h 
and the iLi & Draine (2001) model. Our estimated dust mass is 
also a factor of two higher than found by B06 using the FIR mea- 
surements alone. This is not surprising, since the RJ part of the 
SED, which is sampled at the LABOCA and MAMBO-2 wave- 
lengths, is more sensitive to the mass of dust and less sensitive to 
temperature. Draine et al. (2007) comment that the models used 
in B06 to estimate dust masses are only accurate to within a fac- 
tor of 1 .5 when no submm data are available. Thus, the additional 



submm/mm data points are needed to better trace the bulk of the 
dust mass. 

Adopting a Galactic gas-to-dust ratio (165: 10120051) . our 
dust mass implies a total gas mass of 2.6x10^ Mq. We can com- 
pare this value to the total (molecular atomic) gas mas s us- 
ing measurements from the literature. iBaiaia et all d 19841) find 
an Hi mass of 3.3x10*^ (scaled to our adopted distance of 
9.4 Mpc). Preliminary results from Se ngupta et al. (2003) are in 
good agreement with this value. Bai aia et alj dl991l) used their 
detection of CO at one position in the ring to estimate an up- 
per limit to the molecular gas mass of 4.6x10** (scaled to 
9.4 Mpc), though they note that this value could be uncertain by 
up to a factor of five. Thus, the total gas mass implied by our dust 
mass is consistent with measurements of the total (molecular + 
atomic) gas mass. 



5. Conclusions 

We have observed the Sombrero galaxy (NGC 4594) at 870 \im 
with the LABOCA bolometer camera at the APEX telescope and 
at 1.2 mm with the MAMBO-2 bolometer array at the IRAM 
30m telescope. 

We have detected the dust ring for the first time at submm 
and mm wavelengths. We have constructed a model of the galaxy 
to separate the AGN and dust ring components. The ring radius 
at both 870 (im and 1 .2 mm agrees well with the radius deter- 
mined from optical absorption and atomic gas studies. 

We find that the 70 |^im-1.2 mm SED of the dust ring is well 
fitted by a single grey-body with /3-2 and a dust temperature 
7'd=18.4K. Using this dust temperature we find the ring to have 
a dust mass of 1.6 ±0.2x10^ M©. For a Galactic gas-to-dust ratio, 
our estimated dust mass implies a gas mass that is consistent 
with measurements from the literature. 
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